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New Ideas in Finite Density QCD 
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We introduce a new approach to analyze the phase diagram of QCD at finite chemical potential and temperature, 
based on the definition of a generalized QCD action. Several details of the method will be discussed, with particular 
emphasis on the advantages respect to the imaginary chemical potential approach. 



In 1997, at the end of the previous century, the 
general wisdom on the possibility of understand- 
ing the behavior of matter in extreme conditions 
from first principles was moderately optimistic. 
The determination of the grand canonical parti- 
tion function coefficients of QCD by the Glasgow- 
Illinois group 1 using the Glasgow rewcighting 
method allowed for the detection of a first order 
phase transition at high baryonic density in the 
strong coupling regime. 

One year later, in 1998, we published a paper 
containing very surprising results. Indeed using 
the microcanonical fermion average approach [2] 
and taking the modulus of the determinant of 
the Dirac operator in the integration measure, we 
were able to reproduce with very good agreement 
[Hj the results obtained with the Glasgow method. 
These results were actually surprising since it was 
well known that the phase of the fermion deter- 
minant should play a fundamental role in the dy- 
namics of finite density QCD. 

In the same year we published another paper 
0] in which we did a detailed analysis of the 
grand canonical partition function and detected 
how rounding effects in the standard routine used 
by both groups to calculate the coefficients of the 
grand canonical partition function from the eigen- 
values of the quark propagator matrix, were the 
responsible for the appearance of a phase transi- 
tion. We proposed also in 0] a solution to avoid 
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the rounding problem and the main conclusion 
was that the phase transition observed by both 
groups was missing. In other words, the appar- 
ent progress in the understanding of finite density 
QCD was fictitious. 

From our previous experience we learned that, 
at large fermion masses, the region of the chem- 
ical potential in which the phase of the fermion 
determinant becomes almost irrelevant enlarges, 
and then numerical approaches can be success- 
fully applied. In this way we were able to predict 
a tentative phase diagram for fat QCD 5 which 
confirmed the theoretical expectations. 

This was more or less the general situation of 
three color QCD until two years ago, when an im- 
portant step forward in the field was possible by 
using two alternative schemes. The first scheme, 
based on extracting information on the phase 
structure of the model at finite /i from numeri- 
cal simulations at imaginary chemical potential, 
was implemented in The second scheme [7] 
relies on a two-dimensional generalization of the 
Glasgow reweighting method, and has also been 
combined with a Taylor expansion in the chemical 
potential which allows simulations at large lattice 
sizes in the small /i/T region 

We want to introduce now a new approach 
jH] to simulate QCD at finite temperature and 
baryon density. Even if the theory scheme can 
resemble in some aspects the imaginary chem- 
ical potential approach, its range of applicabil- 
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ity seems much wider. The results of a test of 
our method using as toy model the Gross-Neveu 
model at finite ^t, as well as those for four flavor- 
QCD, can be found in the talk by Angelo Galante 
at this conference. 

Our theory scheme is based on the introduction 
of a generalized QCD action 

S(x,y) = S PG + S T {x,y) + 

1 3 

2 ^Vi(n) (U n>i 1p n+i - U+_ iti 1p n -i) (1) 
n,i=l 

which depends on the two independent parame- 
ters x, y trough 

S T (x,y) = 

-x^i) n r]Q{n) (U nfi ip n+0 - U+_ 0>0 ip n _o) 

n 

+ \y "4>nVo{n) {U n fi^ n+ o + U+_ ofi 1p n - ) (2) 
n 

A natural question now is, what should we ex- 
pect for the phase diagram of this model in the 
x, y plane?. The point x = 1 , y = would cor- 
respond to standard QCD at vanishing chemical 
potential. Now assume we are in the scaling re- 
gion but at a physical temperature T lower than 
the deconfining critical temperature. The point 
x = 1 y = will be in the confined phase. If we 
increase now the inverse gauge coupling /?, the 
physical temperature increases and for (3 values 
large enough the point x = 1 will be in the un- 
confincd phase. This strongly suggests the pres- 
ence of a phase transition point in the y = line 
approaching the x = 1 point in this line by in- 
creasing (3 and eventually crossing it for (3 values 
large enough. This argumentation drives us to 
conjecture the minimal phase diagram shown in 
Fig. 1. The solid line would be a line of phase 
transitions that crosses the y = axis at values of 
x larger than 1 for small (3 values. By increasing 
j3 and keeping fixed the temporal lattice extent 
L t , the critical point on the y — axis moves 
toward x = 1 and eventually crosses it. The dis- 
continuous line in this figure stands for the phys- 
ical line x 2 — y 2 = 1 along which one recovers 
standard QCD. The intersection of the solid line 
with the discontinuous one will give us therefore 




Figure 1. Minimal phase diagram in (x,y) plane. 



the critical chemical potential of QCD at a given 
temperature. 

The numerical analysis of this model for real 
values of x, y is not possible because we find again 
the sign problem. However if we take the y pa- 
rameter as a pure imaginary number, numerical 
simulations are feasible. Assuming that the phase 
transition line of Fig. 1 continues to imaginary 
y = iy, the expected phase diagram in the x,y 
plane will show the qualitative structure of Fig. 
2, where we have incorporated the periodic- 
ity of the model. We have also included in this 
figure the line x 2 + y 2 = 1 (dotted line) which 
contains the only points accessible to numerical 
simulations of QCD at imaginary chemical po- 
tential. As stated in the introduction, one can see 
now how this approach has in principle more po- 
tentialities than the imaginary chemical potential 
approach. Indeed by increasing the inverse gauge 
coupling (3, the phase transition line of Fig. 2 
moves approaching more and more the origin of 
coordinates. In some interval (/3 m ,/?M) the tran- 
sition line intersects the dotted line and then a 
phase transition will appear at imaginary chem- 
ical potential. In such a situation, the physical 
temperature is so high that the system is in an 
unconfincd phase for any real value of the chem- 
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Figure 2. Phase diagram in (x, y) plane. 

ical potential. This is in contrast with the fact 
that, within our approach, simulations can be 
performed at any value of T. 

The practical implementation of our method 
consists therefore in performing numerical sim- 
ulations at imaginary y to determine the phase 
transition line of Fig. 2, and after that, to do 
analytical extensions to real y. This can be done 
mainly in two different ways. In the first scheme 
we fix the inverse gauge coupling and the number 
of temporal lattice points i.e., we fix the physical 
temperature, and determine the equation for the 
critical line of Fig. 1 as an expansion in powers 
of y 2 . The intersection point of this line with the 
physical line will give us the corresponding criti- 
cal value of the chemical potential. In the second 
scheme we fix the parameter x and the gauge cou- 
pling, and look for the intersection of the critical 
line in the (3, y plane, obtained again from analyt- 
ical extension of the critical line in the /3, y plane, 
with the physical line. This scheme corresponds 
to fix the value of /x/T, and the intersection point 
will give us the critical temperature at this value 
of fi/T. 

The main points to remark in the imaginary 
chemical potential approach are: i. the sign prob- 



lem is avoided and standard numerical simula- 
tions can be performed, ii. simulations have to 
be performed at high temperature (large 0) i.e., 
in the deconfined phase for real values of /x, and 
iii. the analytical extension is limited to the range 
it. <- 2L 

T 3 

Our approach shares the first point with the 
imaginary chemical potential case. Furthermore, 
and this is a relevant feature, simulations can be 
performed at any temperature and any value of 
li/T. The price to pay for that is one more pa- 
rameter in the numerical simulations. 

A last suggestion, concerning the Glasgow 
reweighting procedure in two parameters space 
[7] , is the possibility of using the fermion determi- 
nant at the critical point of the y — line instead 
of the determinant at fi = in the integration 
reweighted measure. We suspect that hereby one 
could improve the overlaps. 
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